Scaling exponent of the maximum growth probability in diffusion-limited aggregation 
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An early (and influential) scaling relation in the multifractal theory of Diffusion Limited Ag- 
gregation (DLA) is the Turkevich-Scher conjecture that relates the exponent a m in that character- 
izes the "hottest" region of the harmonic measure and the fractal dimension D of the cluster, i.e. 
D = 1 + otmin. Due to lack of accurate direct measurements of both D and a m in this conjecture 
could never be put to serious test. Using the method of iterated conformal maps D was recently 
determined as D — 1.713 ± 0.003. In this Letter we determine a m i n accurately, with the result 
ctmin = 0.665 ± 0.004. We thus conclude that the Turkevich-Scher conjecture is incorrect for DLA. 



Multifractal measures are normalized distributions ly- 
ing upon fractal sets. As such, they present rich scaling 
properties that have attracted considerable attention in 
the last two decades. In this Letter we address the har- 
monic measure of Diffusion Limited Aggregates [fi| , which 
is the probability measure for a random walker coming 
from infinity to hit the boundary of the fractal cluster. 
This was one of the earliest multifractal measures to be 
studied in the physics literature M, but the elucidation 
of its properties was made difficult by the extreme vari- 
ation of the probability to hit the tips of a DLA versus 
hitting the deep fjords. Thus the understanding of its 
scaling properties has been a long standing issue. These 
scaling properties are conveniently studied using the no- 
tion of generalized dimensions D qi and the associated 
f(a) function || Qj. The simplest definition of the gen- 
eralized dimensions is in terms of a uniform covering of 
the boundary of a DLA cluster with boxes of size £, and 
measuring the probability for a random walker coming 
from infinity to hit a piece of boundary which belongs 
to the i'th box. Denoting this probability by Pi(£), one 
considers ||| 



D„ 



e^o q - 1 log t 



(1) 



where the index i runs over all the boxes that contain 
a piece of the boundary. The limit D = lim q ^ + D q 
is the fractal, or box dimension of the cluster. D\ = 
lim g _ > x+ D q and D2 are the well known information and 
correlation dimensions respectively [|[ ^, It is well 
established by now Q] that the existence of an interesting 
spectrum of values D q is related to the probabilities Pi{t) 
having a spectrum of "singularities" in the sense that 
Pi [£) ~ £ a with a taking on values from a range a m in < 
a < a max . The frequency of observation of a particular 
value of a is determined by the function / (a) where (with 
r(g) = (g-l)D,) 

f(a)=aq(a)-r(q(a)) , ^ = a(q) . (2) 

Of particular interest are the values of the minimal 
and maximal values, a m in and a max , relating to the 
largest and smallest growth probabilities, respectively. 



The maximal value a max was a subject of a long con- 
troversy that was settled only recently (cf. || [| and 
references therein). The issue of ot m i n appears to be one 
of the last of the multifractal properties of DLA that has 
resisted settling. This is the subject of this Letter. 

Consider DLA clusters containing n particles of radius 
a/Ao, and denote the radius of the minimal circle that 
contains the cluster as R n . An incoming random walker 
from infinity has some probability to hit any of the exist- 
ing particles of the cluster. Denote the maximal of these 
probabilities as p m ax- The average of this probabilities 
over many clusters of n particles appears to scale as 



'y/Ao 

Rn 



-a min /D 



(3) 



where for the last step we have used the obvious scaling 
law n ~ (Rn I \ZXo) D . Turkevich and Scher have made 
the plausible assumption that the position of the cluster 
particle associated with p rnax is at the outermost tip of 
the cluster. Thus, a scaling relation can be derived by 
stating that upon adding a new particle to the cluster, 
R n will grow by one unit y/Xo with probability p m ax , or 
will not grow at all with probability 1 — p m ax ■ Then 



dR 
dn 



Ao Pr, 



1/D-l 



(4) 



where again the last step stems from the definition of the 
fractal dimension. Equating the RHS of (||) and (||) we 
get the Turkevich-Scher conjecture 



D = 1 



(5) 



We will show here that this conjecture is incorrect sim- 
ply because the position of maximal probability is not 
at the outermost tip of the DLA cluster. In fact, one 
can introduce in analogy to Eq. (^|) a scaling law for 
the probability to hit the actual tip of the cluster (the 
particle which is furthest away from the origin), i.e. 



(pup) 



A scaling law 
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D = 1 



-a tip /D 



Oitip 



(6) 



(7) 
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is then a tautology. We will show that for DLA au p > 

To achieve accurate estimates of a m i n (and in passing 
of aup) we resort to the method of iterated conformal 
maps that was shown to be extremely useful for dealing 
with DLA and related growth processes. The method 
was amply described before, so we just remind the reader 
that it is based on compositions of fundamental confor- 
mal maps which map the exterior of the unit circle 
to its exterior, except for a little bump at e ie of linear size 
proportional to vA- The composition of these mappings 
is analogous to the aggregation of random walkers in the 
off-lattice DLA model. We shall here use the mapping 
introduced in |hJ which produces two square root sin- 
gularities which wc refer to as the branch cuts, and the 
tip of the bump which we refer to as the micro tip. The 
dynamics is given by 



*(")(«,) =$(«-D (<f> Xni e n (w)) 



(8) 



where maps the exterior of the unit circle to the 
exterior of the cluster of n bumps. The size of the n'th 
bump is controlled by the parameter A„ and in order to 
achieve particles of fixed size we have that, to leading 
order, 



A,, 



An 



|$(n-l)'( e ie„)|2 



(9) 



Using the iterated conformal maps it is very easy to 
keep track of where the maximum growth probability is 
located, and where the outermost tip is, as more particles 
are added. Let us assume that at the (n-l)'th growth 
step the site with the largest probability is located at the 
angle 9 max on the unit circle, i.e. for all 9 
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|$(n-i.)'( e i0 mM ,)| - |$(™-i)'(e i9 )| 



(10) 



When we add a new bump in the n'th growth step the 
position of maximal probability may not change (up to 
reparameterization of the angle 6 max ), or move to the 
new bump. We can easily find the rcparamctcrizcd angle 
and determine the new position from 
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(11) 

If p max ,n is located at 9 n we put 9 max — 9 n in the (n + 
l)'th growth step. Similarly we track the position \z\ max 
on the cluster by finding the value 9u v which assign the 
maximal value of |<J>( n '(e i6l )|. We compute pu p there as 

A direct measurement of a m i n and au p is displayed 
in Fig. [l]. From the direct measurement wc find 
OLmin ~ 0.681 while cttip ~ 0.713. Clearly the latter is 
in agreement with (^) while the former is in disagree- 
ment with (||) (taking as a datum the result of Jl2| , 
D = 1.713 ±0.005). 
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FIG. 1: The average of p m ax (upper, black line) and pu P 
(lower, grey line) versus n. The average is over 20 clusters 
of sizes up to n — 100000. From the values of the slopes we 
estimate of a m in = 0.681 ± .014 and ctu v = 0.713 ± .012. 



The direct measurement, while correct, cannot guar- 
antee that very slow convergence of the power laws as a 
function of n may somehow hide an asymptotic identity 
of a m i n and aup. To remove this worry we adopt now the 
scaling function technique of JlJ] to achieve an accurate 
determination of a m i n . In this approach one acknowl- 
edges that Eq. (J3j) may be realized only asymptotically 
for high values of n. For low and medium values of n, 
(p-max), which is a function of the discrete n and of Ao, is 
in fact a scaling function of a single scaling variable, 
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where we have denoted (3 = a m i n / ' D. The difference with 
Eq. (|) is that f s , p is in general not a linear function of 
its argument, except at exceedingly small values of x, 
when n is very large. In Fig. |2| we demonstrate the 
existence of the scaling function and the excellent data 
collapse achieved using it. In the upper panel we plot 
fs,@(Xo,n) for five values of Ao and 1 < n < 250. In 
the lower panel the same data are collapsed using the 
single scaling variable. We draw the reader's attention 
to the following two observations: (i) the data collapse is 
available immediately, even for the smallest values of n 
Q , and (ii) the scaling function is not linear throughout 
the range explored here. Thus the scaling law (||) is not 
obeyed yet for values of n of the order of a few hundreds. 
The set of parameters j3 and 5 which give the best data 
collapse in the lower panel are /3 — .389 and 6 = .505. 
These parameters are used in the lower panel and give the 
estimate a m i n — 0.666 when assuming that the fractal 
dimension is D = 1.713. 

An even more accurate determination of a m i n is 
achieved next. Taking the data collapse as an evidence 
for the existence of a scaling function, we conclude that 
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FIG. 2: The upper panel shows data of {p m ax} versus n for 
Ao = .8, .9, 1.0, 1.1, 1.2. The lower panel shows the same data 
plotted versus the scaling variable defined in Eq. (|l2|). The 
estimated values of (3 and 8 that lead to the best data col- 
lapse, using least squares, are j3 ma x = .389, 8 ma x = .505. We 
therefore have from D « 1.713 that a m i„ w 0.666. 



for any two pairs of numbers (n, Ao) and (n, Ao) that sat- 
isfy the equation 
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1 



it follows that 



,/n 



f 5,0(^0, n) = /5,/3(Ao,n) 



(13) 



(14) 



These equations offer a calculational procedure. We find 
(Pmax) for a given n and Ao, and then for another value 
n seek the value Ao for which (p ma x) is the same. From 
Eq. ( |l3| ) we then deduce that 

1 logAo-logAo 

Oimin = 7: D ] 7 ] 1~ 1 (15) 

2 log(n + 0) — log(n + 0) 

In Fig. |3| we present the results of such a calculation with 




FIG. 3: Data of a m i n estimated using Eq. ( [l5|) with h = n+1. 
The data is fitted with a cubic polynomial. The polynomial 
intersects the y-axis at the value a m in = 0.665. Using upper 
and lower values of 8, 8 — and 8 = 1 we estimate the 
following bounds on the value of a m in, 0.662 < a m i n < 0.669. 
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FIG. 4: The value of R„ as a function of n (upper curve) and 
the position of the bump belonging to a m in (lower curve) in 
a typical DLA 



n = n+1, and 1 < n < 250. Since S is not known a priori, 
we used the value 5 = 0.505 which was extracted from the 
data collapse in Fig. |[ We checked the sensitivity to 6 by 
bracketing the results with 5 = and 6=1 respectively. 
The data in Fig. || correspond to 6 = 0.505. Fitting the 
data with a cubic polynomial and extrapolating to x — > 
we get the value a m i n ~ 0.665. On the other hand if we 
repeat the procedure using the values of < S < 1 we 
are able to bracket the estimate in the interval 



0.662 < a min < 0.669 



(16) 



We thus conclude the analysis with the estimate a m i n = 
0.665 ±0.004. 

Finally, we explain why the Turkevich-Scher conjecture 
(H) fails. The reason is that the points corresponding to 
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Pmax and pup are not at all the same in typical DLA. In 
Fig. H we present the calculated value of R n , computed 
from the position of largest \z\ on the cluster, compared 
with the position corresponding to the maximal harmonic 
measure. We see that the position of maximal probability 
fluctuates wildly, and the fluctuations do not appear to 
go down with the increase in the cluster size. The loss of 
the conjecture (|^) means that there is no clear connection 



between the spectrum of singularities f(a) and the fractal 
dimension of DLA. As said above, the relation (m) is a 
tautology once the existence of the scaling law ^6|) has 
been established [l| . Since the value of au p has nothing 
to do with the edge of the a spectrum, it appears as hard 
to determine it from first principles as to determine the 
dimension D itself. 
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